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Abstract 

It is shown that the diagrammatic cluster expansion technique for equilibrium aver- 
ages of spin operators may be straightforwardly extended to the calculation of time- 
dependent correlation functions of spin operators. We use this technique to calculate 
exactly the first two non- vanishing moments of the spin-spin and energy-energy cor- 
relation functions of the XY model with arbitrary couplings, in the long- wavelength, 
infinite temperature limit appropriate for spin diffusion. These moments are then 
used to estimate the magnetization and spin-spin energy diffusion coefficients of 
the model using a phenomenological theory of Redfield. Qualitative agreement is 
obtained with recent experiments measuring diffusion of dipolar energy in calcium 
fluoride. 
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1 Introduction 



Experimentally measured quantities in spin systems can often be expressed in 
terms of time-dependent correlation functions of spin operators. [1,2] A well- 
known example is the free induction decay lineshape in solids. [3,4] Another 
example is the rate of spin diffusion, [5] the transport of magnetization or spin- 
spin energy by mutual flips of spin pairs having the same Zeeman splitting. 

The calculation of time-dependent correlation functions can be challenging 
both because of the structure of typical Hamiltonians for spin-spin interac- 
tions and because of the non-trivial commutation properties of spin opera- 
tors. Of particular difficulty is the analysis of correlation functions of more 
than two spin operators. These arise in studying the diffusion of spin-spin 
energy, a problem in which interest has been revived by recent experiments 
that directly observed the diffusion of magnetization and dipolar energy in cal- 
cium fluoride. [6, 7] A phenomenological approach developed by de Gennes[8] 
and Redfield[9] to calculate spin diffusion coefficients based on the knowledge 
of the first few moments of the associated correlation functions agrees well 
with experiments on magnetization diffusion. However, because of the diffi- 
culty of calculating moments for systems with long-range interactions, such as 
calcium fluoride, this approach has not been used to study spin-spin energy 
diffusion in such systems, while magnetization diffusion has only been studied 
to lowest order in perturbation theory in the flip-flop (or XY) term of the 
Hamiltonian. [10] 

In this paper we present a diagrammatic technique for calculating the mo- 
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merits of time-dependent correlation functions, allowing a simplified treat- 
ment of the type of problems mentioned above. This technique extends an 
approach originally developed for calculating static, equilibrium averages in 
spin systems. [11,12,13,14] The extension is based on the cancellation of dis- 
connected diagrams, proved in Appendix B. This cancellation greatly reduces 
the number of diagrams one needs to consider, and constitutes the primary 
advantage of the method. 

The method is illustrated through application to the XY model. This model 
was chosen because it contains the simplest Hamiltonian exhibiting the dy- 
namics of spin diffusion - mutual flips of spin pairs. It is therefore expected 
to qualitatively reproduce the behavior of more complicated systems, such as 
dipolar-coupled spins in high field, for which the dynamics is governed by the 
spin-flip process. It is also useful for comparison to the perturbative limit. [10] 
We calculate the first two non- vanishing moments of the spin-spin and energy- 
energy correlation functions in this model for arbitrary couplings, at infinite 
temperature. The expressions are exact in the long-wavelength limit. From 
these moments, analytic expressions for the diffusion coefficients are obtained. 
Choosing the coupling constants in our calculation to be those of calcium 
fluoride gives numerical values in qualitative agreement with experiments, as 
shown in Table 1. The ratio we find for the diffusion coefficients of magnetiza- 
tion and spin-spin energy is, however, a few times smaller than experimentally 
measured. [7] This may be due to our not having considered the full dipolar in- 
teraction, or to the importance of coherences in the quantum state of the spin 
system, [15] which would not be taken into account by the present approach. 

Besides the moment method, other approaches have yielded spin-spin en- 
ergy diffusion coefficients, such as non-equilibrium statistical mechanics[16] 



3 



and classical simulations. [17] However, the assumptions and approximations 
involved were difficult to justify and gave results which were not in better 
agreement with the recent experiments than those found here. Another recent 
calculation [15] for dipolar interactions was limited to the first two orders of 
perturbation theory in the flip-flop (XY) interaction, and gave similar qualita- 
tive agreement with the experiments. The work presented here should therefore 
complement the previous studies. 



2 Model 



The XY-model for N spins on a rigid lattice is 

N 

H^BvI+ir, (1) 

with Ba = (no sum). The latin indices run over all lattice sites and the 
Jf are spin operators defined by their commutation relations = 5^1?, 

where a, (3, 7 is any cyclic permutation of x, y, z. The if = Ij±ilj are raising 
and lowering operators. The combination of operators, IflJ ', generates mutual 
flips of spin pairs which are responsible for the transport of magnetization and 
spin-spin energy (or heat). The coefficients Bij (i 7^ j) are arbitrary. To make 
contact with dipolar coupled spins, we use 

n di P _ 7 2 fr 3cos 2 % - 1 

a a ~ a r 3 ■ y z ) 

Here 7 is the nuclear gyromagnetic ratio, is the displacement between 
lattice sites i and j, and 9^ is the angle between 17,- and the external magnetic 
field Bq, which is taken to lie along the z-axis. We do not include the Zeeman 
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energy in Eq. (1) as it may be eliminated by a unitary transformation to the 
rotating frame. [1] 

The full Hamiltonian for dipolar coupled spins in a strong magnetic field [1] 
contains an additional term — 2 X)^- B* P J?J?, which we ignore here as dis- 
cussed above. A complementary approach which includes this term but is 
perturbative in Eq. (1) has been discussed earlier. [15] 

The quantities of physical interest are correlation functions of the form 

(S(-k,t)S(k,0)) 

Cs(M)= (S(-k,o)S(k,o)y (3) 

where S(k, t) = J2i e tkZi Si(t) is a spin operator or product of spin operators in 
the Heisenberg representation, and k is the magnitude of the wavevector, which 
points along the magnetic field axis. In the specific cases which we consider 
below, Si is either the local magnetization, Si = —jhli, or spin-spin energy, 
Si = at lattice site i. The angular brackets denote averaging over 

an equilibrium ensemble. For most NMR problems, including spin diffusion, 
it suffices to consider T = oo, so that (• • •) = tr {■ ■ •}/ tr{l}. The extension 
to finite temperature is straightforward and will not be considered here. 

Below we will be interested in the moments of the correlation function, Eq. 
(3). Expanding in Taylor series about t = 0, we obtain 

OO 1 



The even moments M^"' are given by 

M (2n) ( ^ n (l\ 2n (S(k,0)[n,S(-k,0)] 2n ) 

= ( " 1} UJ <S(-M)S(M)> ' (5) 
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where [A, B] n = [A, [A, [...[A, B]...]]], with A appearing n times. The sum in 
Eq. (4) involves only even powers of t because the odd moments are zero. 
These expressions may be derived by expanding S(k,t) = e tHt S(k, 0)e~ lHt by 



Following Redfield,[9] one can obtain an approximate value of the diffusion 
coefficient of S from the first two non- vanishing moments, [18] 



where 0:2 and a>4 are certain phenomenological parameters. We show in the 
next section that each moment is proportional to k 2 at long wavelength, so 
that this expression for Ds is independent of k. Eq. (6) is obtained by matching 
the terms of Eq. (4) to a phenomenological decay function of the form fs(t) = 
9s(t)e~^ Ts , where Ts = (/c 2 -Ds) _1 is the diffusion time. The cutoff function 
gs(t) is different from unity only at times short compared to the spin-spin 
correlation time, T s = h/max.(Bij), which is roughly the time required for a 
single spin flip, and is in principle determined by the microscopic dynamics. 

The exact values of the parameters a<m are related to the manner in which 
the cutoff function gs(t) vanishes at high frequency. For example, a Gaussian 
and step-function cutoff give 



the well-known formula e A Be A = J2^=o tJt [A, B] n and putting the result in 
Eq. (3). 



D s = 



1 1 a 4 (M( 2 )) 3 



(6) 



k 2 r s k 2 \al MW 



(-l)"(2n-2)! 



9siu) 




v ^F2 2 "- 2 (n-l)!' 



(7) 



«2n 



(-1)"2 
7r(2n-l) ' 



g s (u) = B(Ts 
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Both values have the same order of magnitude. Here gs(w) is the Fourier 
transform of gs(t)- Since the shape of the cutoff function is not determined 
within the phenomenological model, Eq. (6) can only be viewed as approxi- 
mate. Nevertheless, this shape is not expected to be drastically different for 
the magnetization and spin-spin energy diffusion coefficients, and therefore 
their ratio can be expected to have a weaker dependence on cutoff. 

3 Calculation of moments 

In this section we calculate the second and fourth moments of magnetization 
and spin-spin energy for the XY model. Since we are interested in the long- 
wavelength behavior, we Taylor expand the correlation function, Eq. (3), in 
k. This gives 



where z^ = zi — Zj, and the terms odd in Zij are zero. The 0(k 4 ) term is 
safely neglected as the correlation {Si(0)Sj(t)) is a rapidly decaying function 
of the distance | — r 3 -|. It depends on products of the spin-spin couplings, 
Bij, which are either short-ranged or, in the case of dipolar coupling, decay 
algebraically on a length scale of a few lattice spacings. We will demonstrate 
this explicitly for each moment. The wavelength, A = 2n/k, is taken to be 
much longer than this decay scale. In the calcium fluoride experiments [6, 7] 
it is at least 10 4 lattice spacings. Expanding the commutator in Eq. (5), we 
obtain 



c s (k,t) 



2 E^(0) 2 > 



+ 0(k 4 ), 
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M, 



(2n) 



2E,(^(0) 2 ) 



2n 



E4E 



m=0 



2n 



m 

V / 



-l) m (W m ^(0)W 2n - m ^(0)), (9) 



for n > 1. Here 



/ \ 

2n 



m 



i 2n ' >! ■ is a binomial coefficient, and we have used 

(2n— m)! ' 



2n 



[W,5,-(0)] an = E 

m=0 



/ \ 

2n 



m 
\ ) 



2n—m 



(10) 



Eq. (9) proves the k 2 dependence mentioned in the last section. 

To calculate the moments for the XY model from Eq. (9), one must evaluate 
averages of the form {TC m Sj(0)Ti. 2n ~ m 5^(0)) . We do this using a diagrammatic 
cluster-expansion technique, [11, 12, 13, 14, 2] extended to Eq. (9) with the help 
of a theorem proved in Appendix B. This technique eliminates the need for 
keeping track of the Kronecker deltas that arise from the contractions of spin 
operators, and allows the identification of the most important contributions to 
Eq. (9) at each n. It is based on an ordered cumulant expansion of spin operator 
averages. For completeness, a brief introduction to ordered cumulants of spin 
operators, also known as semi-invariants, [11,12,13,14,2] is given in Appendix 
A. 
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3.1 Magnetization moments 



Let Si = —jhlf, and consider the expression, 

/ \ 



2n 

:^ = (-i) n+1 E4E 

i,j m=0 



2n 



m 



-l) m (H m I*H 2n - m I*}, 



11) 



in the numerator of Eq. (9). A diagram element is associated to each operator 
in this expression as follows. 





k 



:i2i 



ft = e B kiim 



(13) 



The indices k, I are dummies that are summed over, and in practice can be left 
off of diagrams. One then forms all possible topologically distinct, connected 
diagrams from these elements by joining them end-to-end in all possible ways, 
with the open circles for I z inserted at vertices. The diagram elements are 
numbered based on the order in which they appear in Eq. (11). This order 
must be kept track of because of the non-trivial commutation properties of 
spin operators. For example, the diagrams corresponding to (I'TiHT-CHTi) are 
numbered as follows. 



w 1 02 w 3 "4 "5 06 



Each vertex without a circle is assigned a dummy summation index, and each 
vertex with a circle receives the index corresponding to that circle. The circles 
corresponding to i and j must appear at different vertices, since the zfj factor 
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Fig. 1. Diagram contributing to second moment for magnetization 



in Eq. (11) ensures that i ^ j. To each vertex is assigned an ordered cumulant. 
Each interaction line has an interaction coefficient associated with it that has 
the appropriate indices. E.g. the line k — *— 1 receives a factor of B^. The 
analytic expression corresponding to a given diagram is formed by taking the 
product of all the ordered cumulants and interaction coefficients associated 
with it, and summing over all dummy indices without restriction. The sum 
includes a factor of z\- and the appropriate binomial coefficients appearing in 
Eq. (11). 

Of the total set of possible diagrams, many do not contribute. There are no 
diagrams with free ends, as these represent uncontracted spin operators which 
cause the trace to vanish. Each vertex must have the same number of lines 
leaving as entering, since all ordered cumulants with an unequal number of 
raising and lowering operators are zero. Finally, the disconnected diagrams 
vanish, as shown in Appendix B. 

The only diagram contributing to the second moment is shown in Fig. 1. Its 
contribution to Eq. (11) is 
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^=E4E 

ij m=0 



( \ 



(-l) m X 



=E 



4 x 



1 x 4> + 



=E44 [(«* + -»«- + *» + «* - +»«+ - *») 

-2 (<<+*-))<<- + *>> + <<-*+))<<+-*)>) 
+ («- + *»«+-«» + «+-*»«- + «»)" 



E Z ijBij- 



(14) 



The values of the ordered cumulants are ((H — z)) = | and (( — h 2)) = — |, 
as given in Table A.l. 

The denominator of Eq. (9) may be calculated without diagrams, and we 
obtain 



E(W) 2 ) = j- 



(15) 



Inserting these results into Eq. (9) gives 



M 



(2) 



M 



-k 2 z 2 k B 2 k . 



(16) 



We note that, because of translational invariance, we can drop the summation 
over the dummy index k. 

The diagrams contributing to the fourth moment for magnetization are shown 
in Fig. 2. They are calculated in a similar way to those for the second moment, 
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a) 



b) 





c) 




d) 




Fig. 2. All topologically distinct diagrams containing two circles and four interac- 
tion lines. The diagrams shown here arise in the calculation of the fourth moment 
for magnetization as well as that of the second moment for spin-spin energy. The 
analytic expressions for the diagrams are different in the two cases, however. 



so we omit the details. Table A.l shows that most of the fourth and fifth-order 
cumulants are zero, which enables us to consider only a subset of the orderings 
of the diagram elements. The non-zero cumulants at fourth and fifth order 
correspond to vertices with two ingoing and two outgoing lines, with both 
ingoing lines next to each other in the order (same for the outgoing lines). 
The calculation shows that only the diagrams labelled a), b), and c) in Fig. 2 
contribute. The fourth moment is 



M 



(4) 



M 



-4k 2 



/ j z ik D ik 



E Z ik^ik 



(17) 
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3.2 Energy moments 



The expression in the numerator of Eq. (9) for spin-spin energy, corres 
to Eq. (11), is 



2n 



T 2n = (-i) n+1 E 4 £ 

i,j,k,l m=0 



where $ = Ek,(k^i)^ik, and 



We can rewrite Eq. (18) as 



2n 



—l) {n rtjin ri ik ), 



T 2n = 2(T 2 ( + ) +T, 



(-) 

2n 



2n 



^2n — £ Z ij £ 
m=0 



2n 



T 2n - £ 4 E 

i,j,k,l m=0 



( \ 

2n 



m 
\ J 

( \ 

In 



-l) {n+m+1) (H m H i i t ) H 2n - m H\t ) ), 



where use has been made of the formula 



((A + A*)(B + B*)) = 2Re ((A + A*)B), 



for any operators A and B. 
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We associate the following diagram elements with the operators appearing in 
Eqs. (23) and (24). 



?4 +) — k (26) 



— © (27) 

The diagram element for the full interaction, 7i, is the same as in the last 
section, i.e. Eq. (13). Dummy indices such as k will be left off of the diagrams 
as before. 

The calculation of Eq. (18) is similar to that of Eq. (11). In this case, however, 
the interaction lines due to Tt\^ and "H^ receive an additional factor of |, 
because this factor appears in Eqs. (20) and (21). The final result is multiplied 
by the factor 2 appearing in Eq. (22). 

The diagrams contributing to the second moment for spin-spin energy are 
shown in Fig. 2. These diagrams are exactly the same as the ones arising in the 
calculation of the fourth moment for magnetization. However, their meaning 
is different, as now there are no I" operators, and we deal with a different set 
of ordered cumulants. We note that the diagrams at order (2n) for spin-spin 
energy are always the same as those at order (2n+2) for magnetization. 

One can easily see that the diagrams labelled e), /), and g) in Fig. 2 are zero. 
Associated with each of them is the product of ordered cumulants, ((-I — )) 4 = 
Yq. Because this cumulant factor is the same regardless of the order of diagram 
elements, we can move all the diagrams to the left of the second summation 
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d) 




e) 



g) 



h) 







3 






i) 



Fig. 3. All topologically distinct diagrams containing six interaction lines. Diagrams 
for the fourth energy moment are obtained by placing circles with indices i and j 
at vertices in all distinct ways. 

sign in Eq. (18). For example, diagram e) gives 



ij m=0 



( \ 

2 



4 



'-l) m x 



u m 



15 



^E 



4 x 



" 7 

-* 1 2 

' xE 

m=0 



(2* 



Since the sum over binomial coefficients is zero, we have T 2 (e) = 0. 



By direct calculation, it is also easily found that the diagrams labelled a), 
b), and c) are zero. The only diagram contributing to the second moment for 
spin-spin energy is therefore diagram d) of Fig. 2. Eq. (18) then reads 



^ 2 = 2E4E 

ij m=0 



( \ 



m 



(-l) m x 



(29) 



According to Table A.l, ((H 1 — )) = 0. This restricts the possible orderings 

of the diagram elements, since not all vertices with four lines are allowed. 
Therefore, 



^ 2 = 2E 



fl 3 i { 



(30) 



The product of ordered cumulants is the same for each diagram in Eq. (30). It 

is ((H — )) 2 ((+H )) = (|) 2 (-|) = -|- Multiplying by \ for the two circles, 

we obtain 



1 



T 2 = 2 ( --J ^ ) Y, *lBlB% x [1(3) - 2(2) + 1(3)] 



'.it 



E z ijBikBjk- 

ijk 



(31) 



16 



The denominator of Eq. (9) is given by ((H — )) 2 Bfj = \ Bfj. Inserting 
these results into Eq. (9), we obtain 



H -"T ^Bl ' (32) 



where we have used translational invariance to drop one of the summations. 

The types of diagrams arising in the calculation of the fourth moment are 
shown in Fig. 3. To save space, the distinct topologies are pictured without 
circles. The entire set of diagrams at fourth order is obtained by placing two 
circles at the vertices of the diagrams in Fig. 3 in all possible ways. The result 
is straightforward to calculate, and is 

/if ( 4 ) _ J.2 V- _2 R 2 R 2 ^ k ( B ik B jk + gfcg^fc) 9 2 Eij Z^B^B^B^ 

1V1 H — K z ik n ik n jk ~ ZK v-v D 2 7 K 



2 



A; 2 J2iji z ik (6Bf k B jk B kl B ij B il - \%B ik B 2 - k B k iBijBji + WB^B^B^BuB^ 



J2i Bf k 



The sums over the index k are left off, as usual. The first term in Eq. (33) 
comes from diagrams a) and b). Diagram c) is of the same order of magnitude, 
and gives the second term in this equation. Diagrams d) and e) give rise to the 
third term, and are an order of magnitude smaller for short-ranged or dipolar 
coupling. Diagrams g), h), and i) give the last term in Eq. (33) and are another 
order of magnitude smaller. The general guidelines are that those diagrams 
with the greatest number of lines per pair of vertices are the largest. The ones 
with several pairs of vertices joined by only a single line, such as diagrams 
g), h), and i), are the smallest. There are exceptions to these guidelines (For 
example, diagram /) vanishes, for the same reason as does the corresponding 
diagram at second order.), so care must be taken in their application. As for 



(33) 
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the second moment, the diagrams /) vanish, as do diagrams j) and k). 



4 Numerical results for dipolar-coupled XY model 



The results of numerical evaluation of the moments calculated in the last 
section for = B^ p (see Eq. (2)), are given in Table 1. This corresponds to 
dipolar coupling. We have used values of the gyromagnetic ratio and lattice 
spacing for the fluorines in calcium fluoride of 7 = 2.51 x 10 4 rad s _1 Oe -1 
and a = 2.73 x 1CT 8 cm. Because lattice sums can be evaluated numerically 
only for finite lattice sizes, we used finite size scaling to extract the infinite 
lattice limit. The approach to the infinite lattice value is expected to follow a 
power law. For example, if we approximate the sums by integrals in Eq. (16), 



M$^-k 2 j dhB dip (v) 

a<r<L 



const x / r L dr ( — ) r 2 



= cons tx(I-i). (34) 

We performed a least squares fit to a power law of the quantities in Eqs. 
(16), (17), (32), and (33) as a function of lattice size, for both the [001] and 
[111] orientations of the crystal with respect to the external field. We found it 
sufficient to vary the lattice size between 1 and 81 lattice sites on an edge, in 
increments of 2 lattice sites. This gave agreement with Eq. (34) to better than 
one percent. The numbers in Table 1 are the infinite lattice values extracted 
from the scaling analysis. 

Besides the moments, Table 1 gives the values for the diffusion coefficients for 
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Table 1 

Summary of the results for the dipolar coupled XY model obtained from the moment 
method, with recent experimental values for comparison. 



Moments 


[001] 


[111] 




M$/k 2 (xl(r 7 cm 2 /s 2 ) 


-5.59 


-2.21 




M$/k 2 (xl0 3 cm 2 /s 4 ) 


1.56 


0.130 




M^/k 2 (xHTW/s 2 ) 


-2.80 


-1.08 




M^/k 2 (cm 2 /s 4 ) 


76.2 


28.4 




Results for Gaussian cutoff 


[001] 


[111] 


Dooi/Dm 


D M (xl0- 12 cm 2 /s) 


13.3 


11.4 


1.17 


D H (xl0- 12 cm 2 /s) 


21.2 


8.4 


2.5 


Ratio D n /D M 


1.59 


0.74 




T M (xl0~ 6 s) 


13.4 


35.8 




T n (xl(T 6 s) 


42.8 


43.8 




Experiments (Refs. [6,7]) 


[001] 


[111] 


A)oi/-Dm 


D M Ref. [6] (xl(T 12 cm 2 /s) 


7.1 ±0.5 


5.3 ±0.3 


1.34 ±0.12 


D H Ref. [7] (xl(T 12 cm 2 /s) 


29 ±3 


33 ±4 


0.88 ±0.14 


Ratio Dn/DM 


4.1 ±0.7 


6.2 ± 1.1 
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both Gaussian and step- function cutoff (see Eq. (7)), as well as their ratio. We 
find fair agreement with experiments on calcium fluoride for the magnitudes of 
both diffusion coefficients. For magnetization, our value is slightly high, while 
for spin-spin energy it is slightly low. The ratio D n /D M that we calculate 
is about f.6 for the [001] direction, while in these experiments it is between 

4 and 6. Given the phenomenological nature of the theory we feel this to be 
adequate agreement. For the [111] direction, the results are quite different, 
giving a ratio of diffusion coefficients that is less than one. We cannot account 
for this difference but conjecture that it may be the result of neglecting the 
Ising, or I Z I Z , term from the calculation. 

As an additional check for consistency of this theory we have calculated the 
value of the short time cutoff, Ts, using its relation[10] to the moments of 
the appropriate cutoff function in Eq. (7). As Table 1 shows, T s was found 
to be on the order of 10 - 100 /is for the different cutoff functions and crystal 
orientations that we considered. This is consistent with the assumption that 
Ts is related to the spin-spin correlation time given by the free induction decay. 
The timescale associated with this decay in calcium fluoride is approximately 
20 /is with the external field in the [001] direction and approximately 50 [is 
with the external field in the [111] direction. [4] 

5 Conclusion 

We have shown how the diagrammatic technique for calculating equilibrium 
correlation functions in spin systems may be adapted to the evaluation of 
multi-spin dynamical correlation functions. We used this technique to obtain 
exactly the first two non-vanishing moments of the magnetization and spin- 
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spin energy autocorrelation functions of the XY model at infinite temperature 
and long wavelength. The results were used to estimate the magnetization and 
spin-spin energy diffusion coefficients in the case of dipolar coupling, using a 
phenomenological moment method. We found qualitative agreement with ex- 
periments on calcium fluoride for both diffusion coefficients. The ratio of the 
diffusion coefficient for spin-spin energy to that for magnetization was found to 
be greater than one for the [001] orientation of the external field with respect 
to the crystal axes. However, this is not large enough to accurately account 
for the observations. The orientation dependence of the diffusion coefficients 
was also in qualitative agreement for magnetization, but not for spin-spin 
energy. The lack of any experimentally observed orientation dependence for 
spin-spin energy diffusion leads us to conjecture that some additional, possi- 
bly /c-dependent, decay processes may have been at play in the experiment, 
increasing the observed decay rates. Some artifacts of the coherent time evo- 
lution of the spin system could also have been involved, and would not be 
accounted for in the phenomenological model of irreversible decay that was 
used here. Finally, it is possible that the approximation of dropping the Ising 
(I Z I Z ) term was too drastic. A tractable calculation including this term should 
be possible along the lines presented here. Although we focused here on the 
spin diffusion problem, the generality of the technique should allow for wider 
applicability. 
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A Ordered Cumulants 



Ordered Cumulants are used to simplify the evaluation of averages of spin 
operators summed over the lattice. The non-zero elements of a spin operator 
average contain, in general, several operators with the same lattice site index. 
Since operators with different indices commute, we may rearrange them so that 
all operators with the same index are next to each other, and then factor the 
average into averages over operators at different lattice sites, since traces at dif- 
ferent lattice sites are independent. For example, (if 1*1^ I?) — (it Ik) (If If) 
if i ^ k. Since we only consider a Hamiltonian that is invariant under lat- 
tice translations, the averages in the last expression are independent of index. 
A general spin operator average may be calculated by grouping the opera- 
tors by index in this fashion, in all possible ways, taking care to avoid over- 
counting by not including identical groupings more than once. For example, 

EikAk(i^ik) = ZikAk [6ik(i+ir) + (i - s ik )(i+)(i;)} = j2i At ((i?ir) - (i?)(ir)} 

E ik Ak(I+)(Ik)- Defining the ordered cumulants, ((+-)) = (If If) -(If) (If), 
((+)) = (If ), and ((— )) = (If), f° r an arbitrary index i, we obtain J2ik Ak(lflk) = 

((+-» EiAi + «+»«-» ZikAk- 

Generalizing the above example, we define ordered cumulants, also known 
as semi-invariants, [2,11,12,13,14] iteratively in terms of their factorization in 
cumulants of lower degree. Thus, 
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«+»=</+>, 
«-»=<n, 

((i z i + ))=(i z i + )-((i z ))((i + )), 
«/*/-»=(/«/->-«/*»«/-», 
«/+/-» =(/+/-> -«/+»«/-», 

«/ 2 /+/-» = (ri + r) - «/ 2 / + »«/-» - «/ 2 /-»«/+» 

-((i + n)((n)-((n)((i + ))((n), (a.i) 



and so on. Ordered cumulants are related to spin operator averages in an 
analogous way to the relation of cumulants and averages in probability theory. 
The main difference is that the order of the spin operators within the cumulant 
is important due to their non-trivial commutation relations. Using ordered 
cumulants, it is possible to calculate operator averages without restricting the 
summation indices, as shown in the preceding paragraph. 



A list of ordered cumulants up to degree 5 for spin | and T = oo is given in 
Table A.I. We omit cumulants that differ only by a cyclic permutation of their 
operators. In the limit of infinite temperature in which we are interested, the 
density matrix is proportional to unity, and these cumulants are the same by 
the properties of the trace. We note that this cyclic invariance is not a general 
property at finite temperature. Besides cyclic permutations, cumulants differ- 
ing by any other rearrangements in the order of the operators generally have 
different values even at infinite temperature. Finally, cumulants with unequal 
numbers of raising and lowering operators are zero and are not included in 
the table. 
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Table A.l 

Ordered Cumulants for spin 1/2 and T = oo. We use the shorthand notation + 
for I + , — for I~ , and z for P. Cumulants that do not have the same number of 
raising and lowering operators are zero, and are not included. We also include only 
one of each set of cumulants that differ by a cyclic permutation of its operators. As 
discussed in the text, these cumulants are the same at T = oo. 



«*» =o «+- + -*» =o 



«**» =3 «" + + -*» =0 



((+-)) =\ «+-- + *» =o 



({+-*)) =\ {{++--*)) =-\ 



({- + *)) ="3 «"- + + ^» =\ 



((zzz)) =0 {{ - + - + z)) =0 



({zzzz)) =4 ((+ + + "--)) 



«+-**» =0 «+ + - + —» =\ 
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B Cancellation of disconnected diagrams 



We present a combinatorial proof of the cancellation of disconnected diagrams 
for the time-independent Hamiltonian, Eq. (1). To account for time depen- 
dence or an I Z I Z term, it is possible to proceed by the standard method via 
the interaction picture and ^-matrix expansion. [19] However, the inapplica- 
bility of Wick's theorem for spin operators prohibits the factorization of time- 
ordered products into contractions, and we must eventually use the same type 
of counting argument presented here. 

Eq. (9) contains the term 



£ 

m=0 



n 



m 



.lynQr-nSjHrSi). (B.l) 



There are two types of disconnected diagrams which contribute to this term, 
those in which both Sj and Si appear in the same cumulant, and those in which 
they belong to different cumulants. The latter type of disconnected diagram is 
always zero, because the cyclic permutation symmetry of the trace allows us 
to factor all the cumulants to the left of the summation over m in Eq. (B.l). 

The case where both Sj and Si appear in the same cumulant is slightly more 
involved. Consider the subset of diagrams for which I < n interaction lines 
form the part which is not connected to that containing Sj and Si. The / 
interaction lines can correspond to any of the n TVs appearing in Eq. (B.l), 
whose average may be factored outside the summation over m. Depending on 
which ones we choose to factor out, there will be a different number of Ws to 
the right of the operator Sj. If we choose to leave k < n 7Y's to the right of 
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Sj, we can do this in 



m 



— k j 



n — m 



I — (m — k) 



ways. The sum over m in Eq. 



(B.l) for the set of diagrams with I H's factored out is therefore equal to 



£ 

m=0 



/ \ 

n 



m 

V / 



-i) m (n n - m s j n m s l ) 



n—l 



in 1 ) ^2(n n - l - k s j n h A 2 ) 



k=0 

l+k 

x E(-iy 

m=k 



( \ 



n 



m 



m 



— k j 



n — m 



I — (m — k) 



The product of binomial coefficients in this equation is 



(B.2) 



n 



( \ 



m 



n — m 



m m — k I — (m — k) 

n\m\(n — m)\ 



{n — m)\m\{m — k)\k\(l + k — m)\{n — I — k)\ 

n! 

(m- k)\k\(l + k - m)\(n - I - k)\' 



(B-3) 



The only factors that depend on m are ^ m _ k y^ l+k _ m y = \ 



( \ 

I 

m — k 



. The sum 
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over m in Eq. (B.2) is therefore X^= fc (— 1) 



m 



0. This proves the 



m — k 



V 



/ 



vanishing of disconnected diagrams. 
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